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^ Abstract 

The difference between the flow ellipticities of oppositely charged pions, /\v2 = 'y2[7r~] — '^2[7r^], 
measured recently by STAR collaboration at RHIC shows a linear dependence on the event charge 
asymmetry A± = (7V+ - 7V_)/(7V+ + 7V_): Ai;2(A±) = Ai;2(0) + rA± with a slope r > and a 



O 

Oh 
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fS| non-zero intercept A'U2(0) > of order 10 . We point out two novel mechanisms which could 

^~Z^ explain the non-zero value of the charged elliptic flow /\v2 at zero charge asymmetry A± — 0, 

^ i.e., the non-zero positive intercept Ai;2(0). Both effects are due to the electric fields created by 
O 

^ the colliding ions. These fields have quadrupole asymmetry of the magnitude and the sign needed 



c 



to account for the nonzero intercept A'L'2(0) > in the RHIC data. One of the mechanisms also 



^ involves chiral magnetic effect. This mechanism, although negligible at RHIC energies, may become 

o 

i-H important at LHC energies. 
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I. INTRODUCTION 

The recent analysis [U [2] of STAR event-by-event measurements of the difference be- 
tween the eUiptic flows of positive and negative pions, /S.V2 = t'2[7i"~] — i'2[7i"'''], flnds 
an interesting behavior of this quantity as a function of the event charge asymmetry 
A± = (iV+ - A^_)/(A^+ + iV_). As seen in FIG. |l| A^;2 depends linearly on A±, i.e., 



Av2{A±) = Av2{0) + rA± , 



:i) 



with a slope r > and an intercept A'^2(0) > 0. This observed positive slope agrees with the 



3.75- ■ jt" 



^ t 



STAR Preliminary 



-0.1 -0.05 



A 0.1 



0.05 0.1 

Observed A^ 



Au + Au 200 GeV: 30-40% 
0.15 <p <0.5GeV/c 




-0.04 -0.02 0.02 0.04 



FIG. 1: The STAR data [Ij on the elliptic flow difference between the positive and negative pions 
showing a linear dependence on the charge asymmetry A± = (A^+ — N-)/{N-^ + N-) with a flnite 
positive intercept. 

prediction made in Ref. [3, 4J based on the Chiral Magnetic Wave (CMW) [5, 6j. The CMW 
is a gapless sound-like mode of chiral charges (left-handed or right-handed) propagating 
along the direction of the magnetic field, with a dispersion relation 



U) 



±v^k — iD^k^ + • • • 



(2) 



where the sign, and hence the direction of the propagation with respect to the magnetic 
field, depends on the chirality of the charge fiuctuations. The velocity v^ is given by [5j 



N^eB 
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where a is the inverse susceptibility of the chiral charge, and Dl - the longitudinal diffusion 
constant which depends on the microscopic transport properties of a given system. The 
mechanism described in Ref.[3l S] is based on the observation (also made in Ref.[7j) that 
the initial net electric charge Q is a sum of the net charge carried by left-handed (Ql) and 



right-handed (Qr) chiral carriers (quarks): 

Q = Ql + Qr, (4) 

where on average Qr = Ql. In off-central heavy-ion coUisions at RHIC energies the heavy 
ions can create a transient magnetic field of strength as large as 5 ^ 47Tz^/e ^ lO^^G in the 
direction perpendicular to the reaction plane. The CMW moves chiral charges Ql and Qr 
along the axis of this magnetic field according to Eq. ([2]). Since Ql moves in the direction 
opposite to that of Qr^ the net result will be an electric charge excess in the poles of the 
fireball and a depletion in the central region. See FIG. [2] for a schematic illustration of 




Ql^Qr 



FIG. 2: The Chiral Magnetic Wave (CMW) evolution of the initial charge asymmetry (taken to 
be positive in this example) produces a net electric quadrupole moment responsible for a non-zero 

A'U2 = '?;2[7r~] — 1^2 [tt^] (positive in this this case). The collision plane is horizontal in this sketch. 

the mechanism. Such a charge transport via CMW will result in an electric quadrupole 
moment whose magnitude is naturally proportional to the charge asymmetry. Combined 
with the subsequent radial fiow this leads to the charged elliptic fiow of pions, A'U2, linear 
in A± [l3ll4j. A semi-realistic theoretical simulation performed in Ref.[3l|4] qualitatively and 
semiquantitatively agrees with the STAR data, including a non-trivial dependence on the 
impact parameter (centrality). 

Although the CMW mechanism described above accounts for the nonzero positive slope 
in Fig. [Tj the origin of non-zero positive intercept A'^2(0) has not been so far understood to 
the same extent. An argument based on a simple quark coalescence model in Ref.[8j suggests 
that a positive contribution to the intercept may appear due to initial isospin asymmetry. 
In this paper we shall identify two different mechanisms which are more closely related to 
the CMW effect responsible for the nonzero slope. These mechanisms account for the sign 
as the well as the magnitude of the intercept. 



One way to look at the problem is to note that the nonzero intercept suggests the existence 
of an additional source of electric quadrupole moment, approximately independent of the 
charge asymmetry A± and present even in the neutral plasma, i.e., at A± = 0. We point out 
two sources of such a net positive (out of reaction plane) electric quadrupole moment. Both 
involve electric fields created by the heavy ions which are similar in strength to the magnetic 
fields responsible for the CMW. We describe the two mechanisms in the subsequent sections 
and make an order of magnitude estimate of their effect on the charged elliptic fiow A'^2(0). 

II. SOURCE 1: ELECTRIC FIELD QUADRUPOLE 

A highly relativistic heavy ion carries around itself a Lorentz contracted "pancake" of 
the electromagnetic field perpendicular to its velocity. The magnetic field lines are circular 
around the velocity direction and the heavy ions colliding at a finite impact parameter 
create a net magnetic field in the overlapping region along the direction perpendicular to 
the reaction plane. This field is responsible for the CMW. 




FIG. 3: The net electric field profile from the two heavy-ions (left), which leads to a non-zero 
electric quadrupole moment {right) whose magnitude is independent of the charge asymmetry A± 
for small enough A±. The reaction plane is horizontal. 

On the other hand, the electric field from a single heavy ion is radially directed, perpen- 
dicular to the magnetic field at all points in space-time. The superposition of the electric 
fields from the two colliding ions in the overlapping region is described schematically in 
FIG. IsFleft). The key feature of the eflFect we describe is that this electric field has a 
quadrupole component: the field pointing outward perpendicular to the reaction plane is 
stronger than the field pointing outward in the reaction plane due to partial cancellation of 
the in-plane components. Since the plasma (or pre-equilibrium state, such as glasma, at such 
early times) is conducting, there naturally arise net currents moving charges outward from 



the reaction plane: J ^ E. One easily sees that this charge transport induces a net positive 
electric quadrupole moment along the out-of-plane axis (see FIG. [3]) which, combined with 
subsequent radial collective flow, gives rise to A'U2 > as in the CMW. 

Since the conductivity is non-zero also in a neutral plasma, the eflFect is present even when 
A± = 0. It is also clear from charge conjugation symmetry that the conductivity can only 
depend quadratically on the net charge density. Thus the dependence of the contribution 
of the quadrupole electric fleld eflFect to A'U2 in Eq. dl]) on A± is given by A'U2(0) + 0{A']_). 
The term quadratic in A± is negligible for small enough A± {A± in FIG. [l]is < 0.04). 

We emphasize again that the magnitude of the electric fleld is similar to that of the 
magnetic fleld, which indicates that the order of magnitude of the induced quadrupole 
moment should be comparable to the one induced by the magnetic fleld via the CMW. This 
aspect and other similarities make this mechanism (and also the mechanism described in the 
following section) to be a natural possible explanation of the non-zero intercept in FIG. [!} 
if the nonzero slope is due to the CMW. 

Let us make an order of magnitude estimate of the effect. The elliptic flow difference 
between positively and negatively charged particles is of order A'U2 ^ AQ/N^ot where AQ 
is the total charge (in units of e) transported by this mechanism and A^tot is the total 
multiplicity of charged particles. The value of AQ can be estimated as AQ ^ JAtj^ where 
J is the (units of e) current density induced by the electric fleld, tj is the lifetime of the 
current and A is the typical area of the flreball transverse to current during the lifetime of 
the electric fleld, r^. The area transverse to the current at time te is given by the product 
of the transverse size, which for 30 — 40% central events we take to be of order of the radius 
of the nucleus R ^ 7 fm, and of the longitudinal size, which is of order te'. A = Rte- 

To estimate the current and its lifetime we need to consider two cases depending on 
whether the lifetime te of the electric fleld is shorter or longer than the typical (mean free) 
time between collisions rfree- Although the mechanisms of the charge transfer are different in 
the two cases, the relevant product, J tj^ turns out to be similar. As we shall see below, for 
relatively long-lived electric fleld, te ^ Tfree, the magnitude of the current is determined by 
conductivity and J ^ rpree, while its lifetime is as long as that of the electric fleld tj ^ te- 
On the other hand, for short-lived electric fleld, te <C Tfree, the magnitude of the current is 
determined by the lifetime of the fleld J ^ te^ while its lifetime is limited by the mean free 
time Tj ^ Tfree- In either case the relevant product Jtj ~ r^Tfree- 



Let us now estimate the current J . First consider the case when the hfetime of the electric 
field te is much longer than the time between collisions rpree- This is the hydrodynamic 
regime and thus the current is given by Ohm's law J = aE^ where for the conductivity a 
we could use a lattice QCD result from Ref . [9j , 

a = 0.4er J] ql^{).27eT. (5) 

When Te <C rfree we cannot use Ohm's law. But in this case we could estimate the current 
in the collisionless approximation. The contribution of a each quark or antiquark species / to 
the current is given by Jf = qfTifVf (where f = u^d^ 5, u^d^s) where Vf is the drift velocity 
induced by the electric field by increasing momentum density of this component of the 
plasma by tt/ = te eEqfUf. By Lorentz invariance Vf = TTf/wf^ where Wf is the momentum 
"susceptibility" (given by enthalpy e + p in equilibrium). Putting this together, and taking 
into account that, at A± = 0, quark and antiquark of a given fiavor contribute equally to the 
net current, i.e., Je = Jf + Jj = 2Jf, we find: J = J2f=u4,s ^f = ^cEte J2f=u4,s QFril^/wE, 
where the factor 2 in the last equation counts the equal contributions from quarks and 
antiquarks of a given fiavor F. We can estimate Ue and We using equilibrium momentum 
distribution for massless quarks: ue = 9C(3)T^/7r^ and we = 77r^T^/15. Thus we find 

J ^0.26 tecET^ Y^ ql. (6) 

F=u^d^s 

Comparing to Eq. ^ we can write J ^ 0.65(r^T)a£'. Clearly, when te ^ Tfree, the two 
equations should match, from which we could infer a rough estimate of the mean free time 
Tfree ^ 1.5/T implicit in the lattice result Eq. ([5]). Using this as an estimate of the rfree, we 
can write the result in a form independent of which regime we consider: 



jTj^{).neET^TET,,,,. (7) 



Thus 



The electric field strength can be approximated as eE ^ ^ Ze^ /{4:7rR^) where Z is the 
atomic number of the heavy nucleus, R is its radius, and 7 = v^/(2 GeV) is the Lorentz 
contraction factor which we express in terms of the center of mass energy per colliding nu- 
cleon pair. The life- time of the electric field is roughly given by the longitudinally contracted 

6 



thickness of the coUiding nucleus, which is te ^ 2i?/7, while rfree ^ 1.5/T as inferred from 
the lattice result ^ . Putting this together and using Z ?^ 80 and i? ?^ 7 fm we find 

ZaRT _ 4 / T \ / 10^ \ /200GeV\ 

which has the right order of magnitude to account for the intercept A'U2(0) ^ 3. 10""^ in 
Fig.[l} 

III. SOURCE 2: ELECTRIC QUADRUPOLE VIA CHIRAL MAGNETIC EFFECT 

The other source of a net positive electric quadrupole moment arises through a combi- 
nation of both E and B. Consider a P- and CP-odd quantity E - B. From the electric 
field pattern in FIG. |3]as well as the magnetic field pointing upward perpendicular to the 
reaction plane, one easily recognizes that E - B is positive in the upper half region, while 
it is negative in the lower half region as shown in FIG. [4](left). Via the triangle anomaly 
relation for the axial current J^, 

v^ = ^fE^^)^-^' (10) 

this implies a net positive axial charge created in the upper region, and net negative in 
the lower region. In the presence of the vertical magnetic field, the Chiral Magnetic Effect 
[IOHI2] acting on these axial charges induces an electric charge current (in units of e), 

J = ^ [E"^ '''^B , (11) 

from which it is easily seen that the net eflFect is a development of a positive electric 
quadrupole moment. See FIG. [1] (right) for a schematic explanation. To some extent, 
the mechanism is similar to that of the CMW in Ref . [HJ H] , with a main diflFerence being 
that the mechanism is independent of the initial charge asymmetry A±. 

To estimate the contribution of this mechanism to A'U2, we start from the amount of axial 
charge density created by the non-zero E • B during the lifetime te of the electromagnetic 
field, 

^' ^ ^ ( E ^M ^^^^ ^ 0-1 i^E){eB)TE . (12) 




FIG. 4: The profile of E and B that leads to a creation of axial charges via triangle anomaly (left). 
The Chiral Magnetic Effect (CME) acting on these axial charges induces a net electric quadrupole 
moment (right). 

This corresponds to an axial chemical potential of order /x^ ^ Ja/Xa with the axial charge 
susceptibility xa whose magnitude we can estimate assuming it is similar to the vector 
charge susceptibility x known from lattice QCD [13j to be about x ^ 1.0 T^. The CME in 



Eq. (11) then produces a current 

^ = ^ (E^H ^^(^^) ^ 10"'^"' ieE)ieBYTE, (13) 

which would result in charge separation AQ ^ JAtj with A ^ Rte as before. We could 
also expect that the duration of the chiral magnetic current tj is approximately given by 
the duration of the magnetic field tj ^ tb- Putting this together and estimating the electric 
and magnetic fields again as eE ^ eB ^ ^Za/R'^ one finds for A'^2 ^ AQ/A^tot 

where we used te ^ tb ^ 2i?/7 and i? ?^ 7 fm as before. We see that the magnetically 
induced electric quadrupole effect is negligible at RHIC energies compared to the direct 
electric quadrupole effect in Eq. ([9]). However, the direct effect in Eq. ^ decreases with 



energy ^/s^ while the magnetically induced effect in Eq. (14) is energy independent. 

In addition, the magnetically induced electric quadrupole is sensitive to the lifetime of 
the magnetic field r^. For example, if this time turns out to be much greater than our 
estimate tb = 2i?/7 (due to the conductivity of the medium, as in Ref. [14j), that effect 



may be significantly larger than our estimate in Eq. (14) and could possibly compete with 



the direct electric dipole effect in Eq. ^ at sufficiently high ^/s: 



Comparing this with the estimate of the direct electric quadrupole effect ^ we see that, 
while at top RHIC energy the direct electric quadrupole dominates, the two effects could 
possibly become comparable in magnitude (10~^) at LHC energies due to different ^/s de- 
pendence. 

IV. DISCRETE SYMMETRIES 

We end this note by discussing two discrete symmetries, the charge conjugation C and 
the 180° rotation (or reflection) in the transverse plane, R±^ useful for classifying possible 
mechanisms contributing to the observables: the slope r and the intercept A'U2(0). Under 

A^2(0) ^ -A^2(0) , A^ -^ -A^ , (16) 

therefore r is C-even whereas A'U2(0) is C-odd. 

Since QCD is (7-invariant and the source of (7-violation is in the initial condition, one 
concludes that the physics of A'U2(0) must have its origin in the initial charge asymme- 
try, such as the electromagnetic charge of the heavy-ions (or the isospin asymmetry as in 
Ref. [8j). On the other hand, the slope r could possibly receive contributions from other 
effects unrelated to the charge asymmetry of the initial conditions (as, e.g., in Ref. [15j). 

Under the R± rotation (reflection) in the transverse plane, both r and A'U2(0) are even, 
since V2 is i?^-even. The pattern of the electric field is i?^-even, while the magnetic field 
pattern is i?^-odd. In conjunction with the C-parity discussion above, this tells us that the 
intercept A'U2(0) can be linearly proportional to the electric field, but not the magnetic field. 
Indeed, our first source can be viewed as being linear in the electric field £", whereas the 
second source should be considered as an {EB)B ^ EB'^ effect. 

V. SUMMARY AND CONCLUSION 

We identified two possible sources of the charged elliptic flow Av2{A±) at zero charge 
asymmetry, i.e., A'U2(0). The most straightforward source is the quadrupole pattern of the 
electric fields created by the relativist ic heavy ions in the central overlap region. Although 
these fields last a short time, their magnitude is extremely large due to the well-known 
Lorentz contraction. Our order of magnitude estimates suggest that the contribution of this 
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effect at top RHIC energy is comparable with the observed value of A'z;2(0) and that it should 
be inversely proportional to the collision energy y^, dropping by an order of magnitude at 
LHC energies. 

The simultaneous presence of electric and magnetic fields in the overlap region via the 
chiral magnetic effect provides another source of the quadrupole moment at zero charge 
asymmetry. The magnitude of this effect is less sensitive to the collision energy. It depends, 
however, on the lifetime of the magnetic fields, which could be longer than that of the 
electric fields. We estimate that this effect is negligible at RHIC energies but may become 
comparable with the direct electric quadrupole at LHC. 

We stress that the two mechanisms we point out in this paper give rise to the same 
(positive out of plane) sign of the electric quadrupole moment, which makes the prediction 
for the sign of the effect robust. However, a reliable quantitative prediction of the magnitude 
would clearly require a realistic simulation including spatial variations and fiuctuations of 
the electromagnetic fields [TGHIH]. 
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